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Abstract. We present a first set of results concerning stellar evolution of rotating low-mass stars. Our models 
include fully consistent transport of angular momentum and chemicals due to the combined action of rotation 
induced mixing (according to Maeder & Zahn 1998) and element segregation. The analysis of the effects of local 
variations of molecular weight due to the meridional circulation on the transport of angular momentum and 
chemicals are under the scope of this study. We apply this mechanism to low mass main sequence and subgiant 
stars of population I. 

We show that the so-called ^-currents are of major importance in setting the shape of the rotation profile, specially 
near the core. Furthermore, as shown by Talon & Charbonnel (1998) and Charbonnel & Talon (1999) using models 
without /i-currents, we confirm that rotation-induced mixing in stars braked via magnetic torquing can explain 
the blue side of the Li dip, as well as the low Li abundances observed in subgiants even when /x-currents are taken 
into account. We emphasize that p variations are not to be neglected when treating rotation-induced mixing, and 
that they could be of great importance for latter evolutionary stages. 
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1. Introduction 

Whereas standard stellar evolution models (allowing solely 
for mixing in convective regions) represent a good tool 
for a zeroth order description of stellar observations, they 
cannot account for series of abundance patterns observed 
in various locations of the Hertzsprung-Russell diagram. 

While some abundance anomalies are best explained 
in terms of surface features (as is the case of Am stars 
for example; see Turcotte et al. 2000), others probably 
involve deep mixing (as the He overabundances in O stars; 
see Maeder & Meynet 2000 for details) . Rotation induced 
mixing, in the form of meridional circulation, baroclinic 
and shear instabilities, in conjunction with microscopic 
diffusion has been extensively used in order to reduce the 
discrepancies between standard models and observations, 
and has proved to be quite promising. 

However, following rotation induced mixing in stellar 
evolution codes is challenging and, up to now, many of 
the theoretical results from numerical computations have 
assumed some hypothesis that could be too restrictive. 

Send offprint requests to: Ana Palacios ; apalacio@ast.obs- 
mip.fr 



Here, we wish to further examine the role of horizontal 
chemical inhomogeneities in the evolution of meridional 
circulation. 

As first pointed out by Mestel (1953), meridional circu- 
lation is generated by a departure from spherical symme- 
try, and settles to restore it, but in turn generates itself a 
new dissymmetry: the meridian currents transport mate- 
rial with different mean molecular weights from one place 
to another of the radiative zone, so that a non-spherical 
/^-distribution is set up which retroacts on the circulation, 
diminishing its magnitude. In this early description, as- 
suming that the radiative zone does not suffer any local 
turbulent mixing and that the star remains in solid body 
rotation, this retroaction was actually shown to inhibit to- 
tally the circulation in chemically inhomogeneous regions. 

Kippenhahn (1974) pointed out the importance of tak- 
ing into account the changes in molecular weight, partic- 
ularly because a strong V M may represent a barrier that 
hinders mixing. In the models he proposes to explain the 
Am phenomenon, Vauclair (1977) emphasizes also the im- 
portant role of the /x-barrier and underlines the fact that 
it can considerably delay mixing below the convective en- 
velope of slow rotators. This point, and more precisely the 
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retroaction effect of the /^-currents on meridional circula- 
tion still in the case of Am stars, was rediscussed further 
by Moss (1979). 

This idea of the chocking of meridional circulation by 
mean molecular gradients was pushed further by Zahn 
(1992). In the formalism he develops the vertical velocity 
of meridional circulation, U , can be written as the sum of 
two major terms, one describing the effects of rotation in a 
homogeneous star, and the other accounting for the effects 
of chemical inhomogeneities on the circulation (see § 2). In 
this description, he emphasizes the fact that, due to the 
strong vertical stratification in stars, horizontal motions 
are much favoured compared to vertical ones and that 
strongly anisotropic turbulence is likely to set in (such a 
description is also used by Tassoul & Tassoul 1982). The 
effect of that turbulence is to mix chemicals horizontally, 
thus slowly destroying the horizontal inhomogeneities that 
give rise to /u-currents. 

Vauclair (1999) uses simplified expressions for rota- 
tional mixing and couples them to microscopic diffusion 
in order to examine their combined effect in slowly solid- 
body rotating stars. After the build-up of the horizontal p 
gradients that block the 57-currents as in Mestel's descrip- 
tion, she describes a self-regulated process which could 
freeze both meridional circulation and gravitational set- 
tling. 

In the present paper, we present complete treatment 
of rotation induced mixing, including meridional circula- 
tion, turbulence by shear instabilities and microscopic dif- 
fusion. We will describe the effects of such a combination 
of processes on the transport of both angular momentum 
and chemical species in main sequence and subgiant Pop 
I stars. Concerning the later point, we refer to the papers 
of Talon & Charbonnel (1998) and Charbonnel & Talon 
(1999) (hereafter TC98 and CT99 respectively) who in- 
vestigated these kind of processes under the scope of the 
so-called lithium dip problem, but did not include the ef- 
fects of the variations of the mean molecular weight on 
the mixing. 

We will only consider stars that, according to obser- 
vations in open clusters, lie on the hot side of the lithium 
dip. Indeed, as pointed out by TC98, other transport pro- 
cesses for momentum are known to act in lower mass stars. 
In particular, the rotation induced mixing we implemented 
here fails to reproduce the solar rotation profile as inferred 
from helioseismology (Matias & Zahn, 1997), a conclu- 
sion also reached using other versions of detailed rotation 
induced mixing (Chaboyer, Demarque, & Pinsonneault, 
1995). TC98 suggested that the additional mechanism re- 
sponsible for momentum transport in the Sun could well 
explain the rise of Li abundances on the red side of the Li 
dip, and it is the point of view we adopt here by restricting 
our study to the blue side. 

We intend to apply our description of rotational in- 
duced mixing to abundance anomalies in evolved stars 
(red giants) in a forthcoming study and will therefore set 
here the theoretical stage that we will use in this and in 
papers to come. In § 2, we give a detailed description of 



the physics at work, recalling the basis of the theory pre- 
sented in Zahn (1992) and Maeder & Zahn (1998). We will 
then (§ 3) describe the numerical method as well as the in- 
puts used in the models presented. In § 4 we compare the 
models including mixing to standard evolutionary tracks, 
and analyze in more details the effects of the /x-currents on 
the transport in § 5. We will finally compare our results to 
observations of main sequence and subgiant stars in § 6, 
before giving a summary and conclusions in § 7. 

2. Rotation induced mixing 

Meridional circulation takes place in the radiative inte- 
riors of rotating stars where it is generated by the ther- 
mal imbalance induced by the departure from spherical 
symmetry (Eddington 1925, Vogt 1926). Once established, 
this large scale circulation generates in turn advection of 
momentum and thus favours the development of strong 
horizontal diffusion through shear instability, which is the 
mam source of turbulence^] (Zahn 1992, Talon & Zahn 
1997, Maeder & Zahn 1998). On the other hand, verti- 
cal turbulence is partly inhibited due to the high level of 
stratification in this direction. 

Under these conditions of strong anisotropic turbu- 
lence, one can use the hypothesis of the so-called "shellular 
rotation" (Zahn 1992). In this configuration ~ fi(P) for 
differential rotation is weak on isobars. All quantities de- 
pend solely on pressure and can be split into a mean value 
and its latitudinal perturbation 



/(P,0) = /(P) + /(P)P 2 (cos< 



(1) 



where P2(cos0) is the Legendre polynomial of order two. 
The combined effects of meridional circulation and hori- 
zontal turbulence lead to the vertical transport of chemical 
species and angular momentum. 

The transport of angular momentum obeys an 
advection-diffusion equation treated in Lagrangian coor- 
dinates, namely 



P 



d (r 2 n) 



dt 
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r being the radius, p the density and v v the vertical com- 
ponent of the turbulent viscosity. The vertical component 
of the meridional velocity is given by 



U(r) = 
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where L is the luminosity, M the mass, P the pressure, 
C p the specific heat at constant pressure, T the temper- 
ature, V a d, V and the adiabatic, radiative and mean 
molecular weight gradients respectively (Maeder & Zahn 



1 Baroclinic instabilities might set in as well, but we adopt 
the view that they will be less efficient than the shear. This 
point is discussed by Zahn (1992) and Talon & Zahn (1997). 
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1998). g is the modulus of the effective gravity, defined by 
the hydrostatic equilibrium equation 



g = V$+ ^fi 2 Vs 2 



(4) 



<P = -ip/tp and s = £,t/£, p , where ^, £ p and £ T are the 
coefficients in the general expression of the equation of 
state 



dlnP 
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We have considered non-stationarity, foreseeing papers 
to come where the mixing processes described here will 
be further investigated inside RGB stars. The effects of 
non-stationarity on the circulation, which need to be in- 
cluded as soon as nuclear burning takes place in shells, 
appear through the terms 80/ dt and e grav (Maeder & 
Zahn 1998). 

As shown by Chaboyer & Zahn (1992), the vertical 
transport of chemicals through the combined action of 
vertical advection and strong horizontal diffusion can be 
described as a pure diffusive process. Indeed, splitting the 
concentration Ci of the i th element according to Eq. (|l|), 
the advective transport in the rigorous evolution equation 
for the mean concentration 
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and depend respectively on the rotation profile 
and on the mean molecular weight as 
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where Ht is the temperature scale height, K the thermal 
diffusivity and f e = e/(e + e grav ), with e and e grav being 
respectively the mean nuclear and gravitational energy. 

Xfi and are logarithmic derivatives of the radiative 
conductivity x and the total energy e with respect to jj,, 
while derivatives with respect to T are noted as \t and 
et- 

e m (r) = L(r)/M(r) and p m is the mean density inside 
the considered level surface. 

Dh is the horizontal diffusion coefficient (see Eq. |l^). 9 
and A are the relative variations over an isobar of density 
and mean molecular weight respectively, with 
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can be replaced by a diffusive term, with an effective dif- 
fusion coefficient given by 



(6) 



-Deff = 
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(10) 



The vertical transport of chemicals then obeys a diffu- 
sion equation which, in addition to this macroscopic trans- 
port, also accounts for (vertical) turbulent transport, mi- 
croscopic diffusion (see § 3) and nuclear reactions: 
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Microscopic diffusion appears through the velocity term 
(oc E/diff), and c\ represents the variations of chemical com- 
position due to nuclear reactions. D v is the vertical tur- 
bulent diffusion coefficient, taken equal to the turbulent 
viscosity (Zahn 1992). The total diffusion coefficient Dtot 
for chemicals can be written as the sum of three coeffi- 
cients: 



On 



D 



off 



D v . 



(12) 



A is tightly related to the horizontal inhomogeneities of 
chemicals, and obeys a similar law as the concentration of 
the chemicals under the assumption of shellular rotation, 



<9A , jjd^P _ 6 p ^ 
dt dr r 2 h 



(13) 



(see Chaboyer & Zahn 1992 and Maeder & Zahn 1998 for 
details). Even though this equation admits a stationary 
solution for t r 2 /6Dh, we will follow its time evolu- 
tion since it allows for greater stability in the numerical 
algorithms. 

Because of the low viscosity in stellar interiors, verti- 
cal shear due to differential rotation eventually becomes 
turbulent, and we will assume that this instability domi- 
nates, noting that baroclinic instabilities may also grow. 
However, development of the shear instability is possible 
only under certain conditions which are fulfilled when the 
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flow satisfies to both the Reynolds and the Richardson 
instability criteria. 

The Richardson criterion compares the stabilizing ef- 
fect of the entropy gradient to the amount of energy that 
can be extracted from differential rotation. The instability 
condition is in general written as 



Ri 



N 2 



(dw/dz) 



1 

2 < Ricrit — ^. 



(14) 



This criterion is extremely severe and is seldom satisfied 
in stellar interiors. 

However, as first pointed out by Townsend (1958), 
thermal diffusivity can decrease the stabilizing effect of 
entropy stratification. Maeder (1995) used this character- 
istic to show that in actual conditions, there will always 
exist a small enough length scale that will make the fluid 
unstable. The modified Richardson criterion becomes 



Ri 



N 2 



(du/dzy 



(15) 



where T = v£/6K, v is the turbulent velocity, £ the tur- 
bulent length scale, K is the thermal diffusivity, Nj, is 
the thermal part of the Brunt- Vaisala, frequency and N 2 
is related to mean molecular weight gradients. V can be 
viewed as a turbulent Peclet number, the turbulent vis- 
cosity v v — vt replacing the normal microscopic viscosity. 
The turbulent viscosity is then obtained by setting the 
equality in Eq. (15). Further discussion on these topics 
may be found in Canuto (1998), who suggests that, in 
such a context, the critical Richardson number might ac- 
tually be closer to unity. 

In order for the instability to grow, the associated tur- 
bulent viscosity must also be larger than the microscopic 
viscosity v, as expressed by the Reynolds criterion 



u„ < vRe r 



(16) 



where Re c is the critical Reynolds number (~ 10). 

The estimation of the appropriate formulation for the 
Richardson criterion is critical to give a good estimation 
of the turbulent viscosity, and in turn of the turbulent 
diffusion coefficient. Meynet & Maeder (2000) present a 
critical review of the different expressions proposed so far 
concerning this coefficient. In the present work, we shall 
take the expression derived by Talon & Zahn (1997) for 
the vertical component of the turbulent viscosity 



D„ 



8 Ri cl . it (rdn/dr) 2 
5N 2 /(K + D h )+N 2 /D h ' 



(17) 



which adds to Maeder's prescription by also taking into 
account the weakening effect of the strong horizontal dif- 
fusion on the entropy stratification associated with the 
(vertical) p gradients. The factor | includes, in addition 
to a geometrical factor, the coefficient | found by Maeder 
(1995) when deriving the criterion for shear instabilities 
assuming spherical geometry for the turbulent eddies. 

Dh is the horizontal diffusion coefficient and is such as 
Dh 3> D v , thus leading to strong anisotropic turbulence, 



which is one of the basic hypothesis of the derivation pre- 
sented. This coefficient appears directly in Eq. ([l3|) de- 
scribing the evolution of the variations of A, as well as in 
the vertical and effective diffusivity coefficients. Its evalu- 
ation is the weakest point of the theory first developed by 
Zahn (1992). In the present paper, we use an expression 
linking this coefficient to the advection of angular momen- 
tum 



D h = — 
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Ch is a free parameter of order unity (see Zahn 1992) 
which describes the weakening effect of horizontal turbu- 
lence on the vertical transport of chemicals. It appears 
through Dh in both effective (D e g) and turbulent (D v ) 
diffusion coefficients as can be seen from Eq. ([l0]) and 
(|l7j). As D c g oc Ch and D v tx C'h^ 1 , the net effect of vary- 
ing this parameter on the global diffusion coefficient -Dtot 
might not be very important. 

On the other hand, the magnitude of turbulent dif- 
fusion D v depends on i?i cr it (see Eq. [l?]). Changing the 
Richardson criterion as previously suggested, by using a 
larger critical number i?i cr it , will result in an enhancement 
of D v and in turn of the global diffusion coefficient. 

3. Numerical simulations : inputs 

We present and discuss models computed with the 
Grenoble stellar evolution code, STAREVOL. We included 
the transport of angular momentum and chemical species 
due to the combined action of meridional circulation, tur- 
bulence and microscopic diffusion (in the case of chemi- 
cals), following the formalism of Maeder & Zahn (1998) 
as described in the previous section. 

3.1. Numerical treatment of the transport of angular 
momentum 

The 4 th order differential equation for the transport of an- 
gular momentum (Eq. |2|) is split into four 1 st order equa- 
tions. The system is complemented by Eq. ([l3]) in order 
to follow the feedback of variations of the mean molecular 
weight on transport. 

We use the Newton-Raphson method according to 
Henyey (1964) to solve for angular momentum trans- 
port within the evolutionary code. The first boundary 
conditions impose momentum conservation at convective 
boundaries 
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We complement them by requiring the absence of differ- 
ential rotation at convective boundaries 



Or 



= 



for r = r t , r b 
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We also impose A = there, a natural consequence of the 
fact that convective regions are assumed to be chemically 
homogeneous. r t and rj, correspond to the top (surface) 
and bottom (center) of the radiative zone respectively. J-q 
represents the torque applied at the surface of the star. 
Note that momentum conservation is equivalent to U = 
when there is no convective core. 

In the models presented here, we considered a classical 
perfect gas so that 8 = tp = 1 in Eqs. (||), (5) and (^). We 
also used Ch = 1 and Ri c = 0.25, thus not allowing these 
two parameters to vary. 

The transport of chemicals according to Eq. (|ll]) is 
treated via a finite element method. Boundary conditions 
impose conservation of the integrated mass of each ele- 
ment and can be written as follows 



at 



d_ 

at 
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r I r 2 pdr 

rt 



-r 2 p (UdiftCi) ~ r 2 p (D v + D cS ) ^ 
-Met 



dcj 
Or 

for r = r t 



r pdr 



r 2 p (C/diff c?) + r 2 p (A, + D cS ) ^ 

or 

for r = rt,, 



where M is the mass loss rate. 

We emphasize that rotational induced mixing is 
treated in a self-consistent way and that coupling to the 
evolutionary code ensures feedback of transport mecha- 
nisms on structure at all evolutionary phases. 

3.2. Structure equations 

One of the basic effects of rotation on internal structure is 
to break the spherical symmetry of the star. Rigorously, 
corrections due to this departure from spherical symmetry 
should thus be applied when writing the internal struc- 
ture equations. Meynet & Maeder (1997) developed a for- 
malism suitable in the "shellular rotation" context. They 
showed that rotation has little influence on the hydrostatic 
structure of massive stars, while the effects at the stellar 
surface (in particular on the effective temperature) may 
be important. In the HR diagram, these effects may be 
of the same order of magnitude as those due to mixing. 
In the present work, we made the assumption of negligi- 
ble departures from spherical symmetry, and we did not 
include these corrections (i.e., stellar structure equations 
are solved in a standard way). The effects on surface pa- 
rameters such as L and T c g- discussed here are thus only 
due to rotation-induced mixing. 

3.3. Evolution of the rotational surface velocity 

As already mentioned, we only consider stars hotter than 
6500 K on the ZAMS, associated to the hot side of the 
lithium dip in open clusters. We refer to the work of Gaige 
(1993) on the Hyades F-stars to choose the initial velocity 
of our models and further determine the braking to be ap- 
plied according to the temperature. For further discussion 



on rotation in F-stars, the reader is referred to TC98 and 
CT99 where these aspects are more extensively discussed. 

In the present paper, stars have an initial velocity of 
110 km.s -1 and are assumed to undergo magnetic brak- 
ing while arriving on main sequence (see Tab. ||) as sug- 
gested by the early work of Schatzman (1962). This brak- 
ing was proved to be more efficient for thicker convective 
envelopes, that is to say for lower masses at given metallic- 
ity, a trend which is confirmed by the observations. Their 
rotational velocity is further reduced to only a few krn.s" 1 
when they are on the subgiant branch due to the star's ex- 
pansion. 

The braking law adopted follows the description of 
Kawaler (1988) 

1/2 / m \ -i/2 

(fi < fisat) 



dJ 
df 



-KQ 3 
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R Q 
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= < 



R 

R Q 
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1/2 



M 

M 



-1/2 



(19) 



This formulation corresponds to a field geometry interme- 
diate between a dipolar and a radial field (Kawaler 1988). 
Neglecting the evolution of stellar structure and assum- 
ing solid body rotation, this leads to a Skumanich law in 
t -1 / 2 for the surface velocity. It is widely used in the lit- 
erature (Chaboyer et al. 1995, Krishnamurti et al. 1997, 
Bouvier et al. 1997, Sills & Pinsonneault 2000). The pa- 
rameter f2 sat expresses the fact that magnetic field gener- 
ation saturates at some critical value, as shown by various 
diagnostics of magnetic activity (see Saar 1996 and refer- 
ences therein). This saturation is actually required in or- 
der to retain a sufficient amount of fast rotators in young 
clusters, as originally suggested by Stauffer & Hartmann 
(1987). The constant K is related to the magnitude of the 
magnetic field strength. In particular, in the mass range 
we are studying here, the convection zone shrinks signif- 
icantly from the 1.35M Q to the 1.8M Q model, and this 
parameter, which is normally calibrated on the Sun and 
taken to be a constant in all stars (see Bouvier et al. 1997) 
could very well vary here. 

In studies of stars less massive than ~ I.2M0, there 
is a suggestion that f2 sat varies with mass, and a scaling 
in t~J^ w is invoked (Barnes & Sofia 1996, Krishnamurti et 
al. 1997, Bouvier et al. 1997, Sills & Pinsonneault 2000). 
This is required to explain the fact that K stars rotate 
more rapidly than G stars in the Hyades. However, any 
diagnosis depends on the actual model for momentum dis- 
tribution in the stellar interior. For the masses under the 
scope of this study, braking remains much smaller than 
it is for lower mass stars, as reflected by the fact that at 
the age of the Hyades, even the less massive model is still 
rotating at 30 krn.s" 1 . 

Table |l| presents the rotational characteristics of our 
models, as well as the global convective time-scale as de- 
fined by Kim & Demarque (1996) 

»R 

dr/v conv . (20) 
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Table 1. Rotational parameters of the models 



M 




7~conv 


T e s 


V 


Age 


(M©) 




days 


(K) 


(km.s -1 ) 




1.35 


2.3 x 10 4a 


6.0 


6620 


70.3 


a Per 








6630 


43.6 


Pleiades 








6625 


30.6 


Hyades 


1.4 


4.4 x W 4 " 


4.0 


6770 


97.5 


a Per 








6760 


66.2 


Pleiades 








6745 


46.3 


Hyades 


1.5 


5.1 x 10 4i 


3.4 


7080 


100 


a Per 








7045 


83.0 


Pleiades 








6990 


76.0 


Hyades 


1.8 


7.1 x 10 4U 


0.44 


8280 


110 


a Per 








8020 


94.4 


Pleiades 








7625 


84.1 


Hyades 



The scaling of Kfl^ at is clearly not inversely proportional 
to T CO nv Furthermore, braking for stars more massive than 
~ 1.35M is less efficient than it is in lower mass stars, 
while their convective time-scales are lower. We are thus 
lead to believe it is a mass range for which the value of 
K itself varies, in agreement with our knowledge that it 
is also a transition region for stellar activity. 

Let us note that all our calculations were performed in 
the saturated regime, and that our stars follow the evo- 
lution of the mean velocity for the 3 considered clusters 
(a Per, Pleiades, Hyades). It is thus not possible to disen- 
tangle the effect of a variation of S! sat and a variation of 
K. 

3.4. Treatment of the microscopic diffusion 

We included microscopic diffusion in the form of gravita- 
tional settling as well as that related to thermal gradients, 
using the formulation of Paquette et al. (1986). Atomic 
diffusion for trace elements Li, Be, B must be addressed 
with care in the limited range of effective temperature 
on the left side of the Li dip. Indeed, Richer & Michaud 
(1993) showed that, for main sequence stars with an ef- 
fective temperature in the range ~ 6900 — 7100 K, radia- 
tive forces on Li and Be could lead to increased surface 
abundances, at least if large scale mixing is not efficient 
enough to counter balance them. The recent discovery of 
a slowly rotating Li rich dwarf (Deliyannis et al. 2002) is 
nicely explained in this framework. This had lead TC98 
and CT99 to remove microscopic diffusion on those ele- 
ments, in the corresponding temperature range. Here, we 
adopt a different view point. We compute a first series of 
models including microscopic diffusion (but not radiative 
forces) on all light elements. Then, for stars in the range 
of effective temperature between <~ 6900 — 7100 K, we also 
compute models without microscopic diffusion on Li, Be 
and B. The differences between these two series of models 
give an indication of the amount of microscopic diffusion 
that remains in our rotating models. As we shall show, 
this amount is minute. One must keep in mind though 



that the real effect of microscopic diffusion in these stars 
will be to slightly enhance Li and Be abundances rather 
than diminish them. 

Let us note that radiative forces have a magnitude which 
is similar to that of settling (see Richer et al. 2000 Figs. 

5 & 6). Thus, if mixing is efficient enough to (almost) in- 
hibit settling, it will also inhibit levitation by radiation. 
For main sequence stars with T e g > 7100 K, radiative 
forces do not play a major role compared to atomic diffu- 
sion, which is thus the only process (apart from rotation- 
induced mixing) which is applied. 

3.5. Input physics 

A rather detailed description of our stellar evolution code 
can be found in Siess et al. (1997) (see also Sicss et al. 
(2000) for more recent aspects). As far as main sequence 
and subgiant models are concerned, let us briefly sum- 
marize here the relevant physics. Our equation of state 
has been developed from the Pols et al. (1995) formal- 
ism. Thermodynamical features of each plasma compo- 
nent (ions, electrons and photons, as well as H~ and H 2 ) 
are obtained by minimizing the Helmholtz free energy that 
includes separately non-ideal effects (Coulomb shielding 
and pressure ionisation). In input, instead of the usual 
pressure or density variable, it uses a new independent 
variable that describes the electron degeneracy. This al- 
lows in turn to treat analytically the ionisation process and 
provides very smooth profiles of thermodynamical quanti- 
ties. It is particularly well suited to treat partially ionised 
and partially degenerated regions, pressure ionisation and 
Coulomb interactions in a stellar evolution code. More de- 
tails about this equation of state are given in Siess et al. 
(2000). 

Radiative opacities are interpolated from Alexander 

6 Fergusson (1994) at temperatures below 8000 K, and 
from Iglcsias & Rogers (1996) at higher temperatures. 
Atmospheres are integrated in the gray, plan-parallel and 
Eddington approximations. Our nuclear reaction network 
allows to follow the abundance evolution of 53 species 
(from 1 H to 37 C1) through 180 reactions. All nuclear re- 
action rates have been updated with the NACRE com- 
pilation (Angulo et al., 1999). The value of a p — l/H p , 
the ratio of the mixing length to the pressure scale height, 
is taken equal to 1.75. No overshooting is considered for 
convection. 

4. Global characteristics of the models 

We present evolutionary results for five different masses 
(1.35, 1.4, 1.5, 1.8 and 2.2 M Q ) at Z = 0.02. Models include 
complete treatment of angular momentum and chemicals 
transport, taking into account the effects of changes of /x 
{E^ ^ 0) and applying microscopic diffusion to all ele- 
ments. 

These models sample the effective temperature domain 
covered by the observations of main-sequence stars on 
the hot side of the lithium dip in young open clusters 
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Table 2. Characteristics of the stellar models and rotation 
velocities. Models computed with Z = 0.02. The turn-off 
(TO) is defined as the point for which the central hydrogen 
mass fraction is lower than 10 -9 . 
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with nearly solar metallicity such as the Hyades, Coma 
Berenices and Praesepe. They also span the mass range of 
Pop I subgiants to which our results will be compared in 
a further section (§ 6). 

In Fig. [j], we plot the evolutionary tracks of our ro- 
tating models in the Hertzsprung-Russell diagram (solid 
lines) , comparing them to their standard analogs (dashed 
lines) . 

Table |^ presents some characteristics of the computed 
models with and without rotation. All models were com- 
puted with the same set of parameters for rotation induced 
mixing (C/, = 1 and Ri cri t = 0.25), but with different ro- 
tational histories (see § 3.2). 

One can see clearly how mixing affects main-sequence 
evolution. In rotating models, the turn-off is moved to- 
wards cooler temperatures while the luminosity is slightly 
enhanced. Due to surface effects, still cooler effective tem- 
peratures (by up to ~ 100 K in the fastest stars considered 
here) would be obtained if the effects of rotation on struc- 
ture equations were taken into account (Perez Hernandez 
et al. 1999). Let us note that the actual value depends on 
the inclinaison angle of the star. 

Our models are computed without convective core 
overshooting. We checked however that for the 1.8M Q 
model, the effect of mixing on main sequence lifetime and 
on the turnoff in the HRD is equivalent to an overshoot- 
ing distance for convection of about 0.1-ffp (where Hp is 
the pressure scale height at the edge of the classical core). 
This value is smaller than the one (0.2) typically used to fit 
the observed color-magnitude diagrams of young and in- 
termediate age galactic open clusters (Meynet et al. 1993, 
Lebreton 2000 and references therein) . Whether the ro- 
tating models would fit the cluster isochrones remains to 
be checked. This requires a more extended grid of masses 
and is out of the scope of the present paper. 

On the subgiant branch, stars evolve at higher lumi- 
nosities. This is due to the effect of microscopic diffusion 
during the main sequence phase. Indeed, even if it is partly 
balanced by the action of large scale mixing, microscopic 




3.6 



l°g(T„,) 



Fig. 1. HR diagrams for the five masses considered in this 
paper. Dashed lines are for standard models and solid lines 
are for models with the "complete mixing" and Ch = 1. 
Masses are as indicated on the plot. Z = 0.02 for all models 
presented. 

diffusion leads to the build up of a negative helium gradi- 
ent at the base of the convective envelope. When the first 
dredge-up begins, the envelope encounters regions with 
a lower opacity (due to enhancement of helium) leading 
to enhanced luminosities compared to standard models. 
When reaching the RGB, differences in effective tempera- 
ture and luminosity between standard and rotating models 
become marginal: both tracks are almost identical in the 
HR diagram (see Fig. |l|). 

Models with rotational mixing behave as if their metal- 
licity was lower than in the standard cases. This trend was 
already emphasized by CT99. 

On the other hand, the lifetimes of rotating models 
are enhanced with respect to those of standard models. 
In fact, due to the positive slope of the hydrogen profile, 
mixing slightly feeds the core with fresh hydrogen fuel. 
The exhaustion of hydrogen in the central region is delayed 
and the time spent on main sequence increases. 

In the mass range we are studying here, the less mas- 
sive the star, the more turn-off ages differ between rotat- 
ing and non-rotating models^ This trend with mass has to 
be related to braking, which is more important for lower 
mass stars as imposed by both theory and observations. 
Stronger mixing corresponds to a more efficient extraction 
of angular momentum at the surface, thus leading to en- 
hanced meridional circulation and mixing. There should 

2 The time spent on the main sequence is increased by more 
than 20% in the lower mass stars, and by ~ 10% in the higher 
mass models. 
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thus be a correlation between mixing and braking, as we 
will show from the comparison of our models with obser- 
vations (see § 6). 

5. Effects of V M and ^-currents 

To discuss the effects of chemical composition variations 
on transport, we will focus on the 1.5 Mq model. We 
present a comparative study of the results obtained with 
and without /i-currents (E M = or ^ 0. In both cases 
all parameters are kept identical in order to make the com- 
parison clearer. Furthermore, we applied a torque so that 
surface velocity at the age of the Hyades (~ 700 Myrs) is 
the same whether is zero or not. 

5.1. Meridional circulation velocity U(r) 

Figure || presents profiles of U, £1 and the diffusion coef- 
ficient Z?tot at different times on the main sequence, in- 
dicated in terms of central hydrogen mass fraction X c . 
We superimpose the profiles obtained for E^ = (dashed 
lines) and E^ ^ (solid lines). 

When neglecting /i-currents (E^ = 0), the meridional 
velocity is negative in the entire mixing zone. There is only 
one meridian loop that brings matter upwards at the equa- 
tor and down in polar regions. This leads to the extraction 
of angular momentum in the radiative zone, a consequence 
of braking. The amplitude of meridional velocity decreases 
as the star evolves and braking slows down. 

When E^ is not zero, some changes appear. There are 
now two loops of circulation. In the outer parts of the 
radiative zone, the negative loop keeps on transporting 
material from the equator up along the pole, in response 
to the extraction of angular momentum. In the interior 
however, circulation is positive, revealing an inward trans- 
port of angular momentum. Equations (^), (5) and (Q) 
already suggest a strong interaction between £1- and /it- 
currents. Apart from the first term in Eq. (|^), which in- 
cludes Eddington-Sweet and Gratton-Opik terms, Eq and 
E^ are almost mirrored. At the beginning of the evolu- 
tion, \Eq\ > \Efj,\, for /i-currents have not settled yet. As 
the star evolves and //-gradients built-up, //-currents may 
grow. \Efj,\ and \Eq\ eventually become of the same order 
of magnitude, and U disminishes. 

Contrary to what was expected by Mestel (1953) the 
circulation is not frozen. Indeed, the two terms E^ and 
Eq, while being of the same order, never compensate each 
other as can be seen in Fig. |^, where U becomes very 
small, but never disappears completely. Indeed, due to 
the strength of the horizontal diffusion, the circulation 
must remain to continuously refresh the horizontal // fluc- 
tuation. However, the difference between both contribu- 
tions Eq and E^ decreases significantly as the star evolves. 
Meridional circulation velocity is thus quite slow near the 
end of the main sequence, but transport of chemicals still 
occurs, as can be seen from the diffusion coefficient in 

Fig. |. 



5.2. Rotation profile S7 

The inner transport of angular momentum also affects the 
rotation profiles (Fig. |^ middle panels): when including //- 
currents, differential rotation grows larger. In particular, 
the core gains momentum and gets accelerated as evolu- 
tion proceeds instead of keeping a nearly constant rotation 
rate. This can be of crucial importance when further fol- 
lowing the evolution of the star on the red giant branch. 
This behaviour may be interpreted as follows: the rate of 
differential rotation and the relative variations of mean 
molecular weight are related to each other via the equa- 
tion of state (Eq. ||) . Observations and theory show us that 
the circulation in main sequence stars on the hot side of 
the Li dip must act to extract angular momentum. In this 
configuration, the term Eq creates an anti-clockwise cir- 
culation transporting matter from the pole down to the 
equator. E^, on the otherside, opposing Eq (see Eq. ||), 
generates a clockwise circulation. The net effect is a re- 
duction of angular momentum transfer to the surface and 
an enhancement of differential rotation. 

In Fig. [| we present profiles of 9 and A at given evo- 
lutionary points. The bold line represents the 9 profile in 
the absence of //-currents. The difference between both is 
striking. It shows, as was already pointed out by Talon et 
al. (1997) in their study of a 9 M Q star, that the descrip- 
tion of the rotation profile and of the circulation is not 
anymore a matter of 9 (or fl) alone, but of 9 — A. Indeed, 
9 is then conditioned by A due again to their mirrored 
evolution (Eqs. 5 and 0). 

The large central differential rotation is a key ingredi- 
ent in core evolution. Indeed, //-gradients in that region 
inhibit the development of a large turbulent diffusivity 
(see Eq. [l3]). However, the large shear resulting from the 
mirror evolution of 9 and A significantly increases local 
core mixing^]. 

Let us mention the work of Theado & Vauclair (2001), 
who performed calculations including the effects of fi- 
currents according to Vauclair (1999) in order to re- 
produce the low dispersion found in the "Spite lithium 
plateau". Assuming slow, constant, solid body rotation (a 
few km.s" 1 in typical halo dwarfs), and that the transport 
of chemicals is progressively blocked when E^ = Eq (a sit- 
uation they call creeping paralysis), they achieve a disper- 
sion of 0.1 dex in the lithium plateau, in agreement with 
observations. We point out the fact that their scenario 
implies solid body rotation and thus prevents the rotation 
profile to react to horizontal molecular weight gradients as 
is the case in our description. Solid body rotation would 
be enforced in the presence of a small magnetic field or 
by gravity waves (Talon & Charbonnel, in preparation), 

3 As explained in Talon et al. 1997, in a non-evolving model, 
Eq. (2) has a stationnary solution in which advection of mo- 
mentum by the circulation is compensated by turbulent trans- 
port. There is thus a O equilibrium profile. When adding fi- 
currents, about the same equilibrium exists for O — A. This 
leads to following the behavior of A, resulting in large dif- 
ferential rotation in the core. 
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Fig. 2. Profiles of the vertical component of the meridional velocity U r (left panels), angular velocity £1 (middle panels) 
and diffusion coefficients (right panels) inside a 1.5 M© , Z = 0.02 star. Solid lines are for the model with ^ and 
dashed lines for the model with = 0. Central hydrogen content is as indicated. In the left panels, bold parts in the 
solid lines indicate positive values of U r . 



which might be the case in halo stars as well as in the 
stars lying on the cool side of the Li dip. However, as dis- 
cussed by TC98, such an additional process of extraction 
of angular momentum is not expected to be efficient in 
the more massive stars lying on the left side of the Li dip 
which are under the scope of this study. 

5.3. Dif Fusion coefficients 

The examination of the total diffusion coefficients (from 
Fig. ||), reveals that, except towards the centre, they ap- 
pear to be very much the same, whether E^ is zero or 
not. The transport of chemicals, contrary to that of angu- 
lar momentum, will thus not be significantly affected dur- 
ing main sequence evolution (see next section). However, 
while the total coefficients are of similar amplitude, they 
are not the result of the same combination of effects. Fig. ^| 
presents various diffusion coefficients in both cases for a 
central hydrogen content X c = 0.2. We can firstly empha- 
size the fact that the shellular rotation is actually realistic, 



for Dh 3> -D ff and D v . Secondly, when ^ 0, the tur- 
bulent diffusion coefficient D v is larger than the effective 
coefficient D e g almost everywhere in the radiative zone. 
This is a consequence of the enhancement of the differen- 
tial rotation rate in this case, which favours the growth 
of turbulence^. Due to the strong anisotropy of turbulent 
motions, the buoyancy restoring forces in the vertical di- 
rection are slightly reduced because of the horizontal diffu- 
sion, and vertical turbulent mixing is allowed. The profile 
of the effective diffusion coefficient is influenced by U (r) 
(Eq. |l(^) , so that it reflects the loops of the circulation in 
the presence of ^(-currents, each "drop" corresponding to 
an inversion of the direction of the flow. 

When neglecting variations of /i, turbulence can only 
marginally develop in a small portion of the radiative 
zone and the vertical diffusion coefficient D v is reduced 
to molecular viscosity which is always smaller than D e g. 



4 This appears clearly when turning to Eq. (p"7|), through 
which D v oc (dQ/dr) 2 . 
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Fig. 3. Vertical component of the meridional circulation velocity, U as a function of the reduced radius at the same 
epochs inside the same star. 
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Fig. 4. Profiles of 8 and A inside the 1.5 M Q star for three 
different values of the central hydrogen content. Solid lines 
are for 8, dashed ones for A. Thin lines stand for models 
where "complete mixing" (E M ^ 0) is applied while the 
bold line gives the profile of 8 when E M = 0. 

In the central part, we may notice that the total dif- 
fusion coefficient is lowered by one order of magnitude 
in the presence of /x-currents. As pointed out before, /> 
gradients in that region inhibit the development of an ef- 
ficient turbulent diffusivity. This appears clearly in Fig. 



(left panel), where we see D e g dropping down below 10 1 
just above the core. D v also decreases in the same re- 
gion, where contrary to the rest of the radiative zone, it is 
dominated by the molecular viscosity rather than by the 
turbulent viscosity (there is no turbulence). 

Figure |^ presents a comparison between the evolution 
of the lithium surface abundance in the 1.5M Q rotating 
star when considering or not the effects of /x-currents. Both 
tracks are similar, but the depletion of lithium during the 
main sequence is larger in case of — 0. The surface 
abundance of lithium differs by more than 0.3 dex at 700 
Myrs, which is more than the observational errorbar com- 
monly assumed. Accounting for \x variations reduces the 
effect of mixing when all others assumptions remain the 
same. We have seen that /i-currents rule the evolution of 
the rotation profile. They lead to an enhancement of the 
degree of differential rotation, which in turn allows for the 
horizontal turbulence to develop more easily. This hori- 
zontal turbulence hinders the transport towards the ver- 
tical direction. This effect shows up through the surface 
abundance of lithium, that is less depleted when ^ 0. 
Nevertheless, this abundance feature can not be advocated 
as an observational clue to the action of the /i-currents. 
Indeed, braking plays a major role in determining the am- 
plitude of mixing (and thus the amplitude of lithium de- 
pletion), and thus may rub off the effects of the variations 
of \i. The surface velocity reached for T e g ~ 4900 K is of 
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Fig. 5. Profiles of the different components of the diffusion 
coefficient and of the thermal diffusivity as a function of 
the reduced radius when ^ (left) and — (right). 
The models correspond to the 1.5M Q star with X c = 0.2. 



Fig. 6. Evolution of the surface abundance of lithium with 
temperature in the 1.5 star for the two different ap- 
proaches of mixing as mentioned on the plot and in the 
text. 



a few km.s -1 in both cases, velocity for the model with 
Efj, = being slightly larger (9 km.s -1 versus 6.5 krn.s" 1 
when E^ ^ 0), because of the different braking constant 
used in this case in order to achieve the same surface ve- 
locity at 700 Myrs as in absence of /j,-currents. Both veloc- 
ities reached at the end of the computation are consistent 
with observations (Gaige 1993). As a result, accounting 
for the effects of /i-currents under the assumptions de- 
scribed in § 3 does not make any significant change in 
the outcoming surface abundances compared to the case 
without variations of /i. The observed abundances used as 
constraints do not allow to distinguish between these two 
cases. However we pointed out differences in the final rota- 
tion profiles, which might change significantly rotational 
transport during later evolutionary phases. 

6. Comparison with observations 

We now compare the predictions of our rotating-models 
with the observations of Li, Be and CNO on the hot side 
of the so-called Li dip at the age of the Hyades, and of Li 
in subgiants of solar metallicity. 

6.1. Li, Be and CNO on the hot side of the Li dip in 
the open clusters 

The drop-off in the Li content of main-sequence F stars 
in a narrow range of effective temperature centred around 
6700 K, or "Li dip" (Wallerstein et al. 1965, Boesgaard 
& Tripicco 1986, Balachandran 1995), is the most striking 
signature of transport processes that occur inside these 



stars. We refer to TC98 and CT99 for a detailed discus- 
sion of the different processes which have been proposed to 
account for this feature. These papers show how rotation 
induced mixing shapes in a very natural way the hot side 
of the Li dip. As already mentioned, we did not consider 
the effects of the [i gradients in these previous papers. 
However, as discussed in § 4.2.3, the total diffusion coeffi- 
cients, and thus the resulting transport of chemicals, is not 
significantly affected by these terms during main sequence. 
As a consequence, the comparison of the predictions of Li 
and Be surface abundances with observations on the hot 
side of the Li dip remains highly satisfactory, as can be 
seen in Figs. |j| and ||. 

In CT99, microscopic diffusion on LiBeB had been 
turned off (while it was applied to all other chemicals) for 
stars with T e g higher than 6900 K, because of the competi- 
tion between radiative acceleration and gravity (Richer & 
Michaud 1993, Richer et al. 2000, see discussion § 3.4). In 
Figs. ^ and || we present the results of computations with 
and without gravitational settling of LiBeB at the effective 
temperature where radiative forces can over-compensate 
gravitation. As can be seen, for the rotationnal velocities 
considered here, the differences are minute. In addition, 
the effect on stellar structure itself is null, since LiBeB 
are trace elements. 

Observations of C, N and O are precious clues on trans- 
port mechanisms inside main-sequence F-type stars, and 
on the competition between atomic diffusion and macro- 
scopic processes. Indeed, in the case of pure atomic diffu- 
sion, these elements are expected to be under abundant in 
the Li dip, as shown in Fig. 0. On the other hand, in our 
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Fig. 7. Comparison of the models (hexagons and open 
bold star) with observations in the Hyades (crosses), 
in Praesepe (squares) and in Coma Berenices (circles 
and triangles). Downward triangles are for upper limits. 
Observations are from Burkhart & Coupry (1998,2000) 
and Boesgaard (1987). The corresponding masses of the 
models are as indicated on the figure. All models presented 
have Z = 0.02 and an initial velocity of 110 km.s^ 1 ; the 
reader is referred to Tab. || for details on braking. For the 
1.5 Mq model, open bold star is obtained when not in- 
cluding microscopic diffusion on LiBeB while for hexagons 
microscopic diffusion was included also for these light el- 
ements (see text for further details). 



models rotation induced mixing strongly inhibits the set- 
tling effects. As a result, the surface abundances of CNO 
do not vary in the T e g range considered. This is in per- 
fect agreement with observations in the Hyades (Varenne 
& Monier 1998, Takeda 1998; see also Boesgaard 1989, 
Friel & Boesgaard 1990, Garcia Lopez et al. 1993), as can 
be seen in Fig. |[ Let us stress that in the temperature 
domain of A-type stars, most of the observations report 
abundances for chemically peculiar objects, namely Am 
stars, that are slow rotators in which the action of gravi- 
tational settling is important on CNO elements. These pe- 
culiar patterns might be achieved by considering the com- 
bined effects of large scale circulation, turbulence, gravita- 
tional settling and radiative forces (the effects of the latter 
not being included in our evolutionary code) in models 
with a slow initial rotation. Such computations are not 
presented here for this is out of the scope of the present 
study to focus on such a particular phenomenon. 



Fig. 8. Beryllium abundance as a function of effective 
temperature in the Hyades (Boesgaard & King 2002). The 
pentagon, hexagon and triangle are the values obtained for 
the 1.35 M©, 1.4 M and 1.5 M Q models respectively. The 
empty triangle is a 1.5 Mq model where microscopic dif- 
fusion is not applied to LiBeB. Refer to Tab. 1 for further 
details on braking. 

6.2. Li in subgiants 

Early observations by Alschuler (1975) of a few field stars 
crossing for the first time the Hertzsprung gap indicated 
that lithium depletion starts earlier than predicted by 
standard dilution for stars more massive than 2 Mq . Later 
on, significant lithium depletion was observed in a non 
negligible number of slightly evolved field stars (Brown et 
al. 1989, Balachandran 1990). do Nascimento et al. (2000) 
studied the behaviour of lithium in a sample of field sub- 
giants observed by Lebre et al. (1999) for which Hipparcos 
data allowed the precise determination of both evolution- 
ary status and mass. They confirmed that stars originating 
from the hot side of the dip present a large range of lithium 
abundances which can not be explained by standard dilu- 
tion alone and which reflect different degrees of depletion 
of this element while on the main sequence, even if its sig- 
nature does not appear at the stellar surface at the age 
of the Hyades (see Vauclair 1991, Charbonnel & Vauclair 
1992). In this cluster, dilution is not sufficient to explain 
lithium values in evolved stars (which have masses of the 
order of 2.2 M Q ), while main sequence stars present galac- 
tic abundance (Boesgaard et al. 1977, Duncan et al. 1998). 
This has been observed also in open clusters with turnoff 
masses higher than ~ 1.6 M Q (Gilroy 1989, Pasquini et 
al. 2001). 

In Fig. ^ we compare our predictions for the surface 
Li abundance of various models with the observations in 



Palacios et al.: Rotational mixing in Pop I low-mass stars 



13 



i o 
u 

-0.5 
-1 
0.5 

sT o 
z 

-0.5 
-1 
0.5 
X 

-0.5 
-1 



• * # 



7r TT 



- 



X 



I 



H 1 h 



I 



7000 
T 

1 off 



z 



Fig. 9. Abundances of carbon, nitrogen and oxygen vs ef- 
fective temperature in the Hyades, presented as asterisks 
and open stars for A and Am stars respectively (the lat- 
ter ones are not meant to be reproduced here; see text 
for details). Observations of C and O are from Varenne 
& Monier (1998). Observations of N are from Takeda et 
al. (1998), crosses are derived from the 8629 A line and 
stars are from 8683 A line. Abundances are given rela- 
tive to solar values from Grevesse & Noels (1993). The 
typical errorbar for observational values is given for each 
panel. Models with complete mixing, a turbulence param- 
eter Ch — 1 and braking as reported in Tab. || are repre- 
sented by the filled points with hexagon, pentagon, square 
and triangle standing for 1.8, 1.5, 1.4 and 1.35 M Q respec- 
tively, with Z* = 0.02. Lines are from calculations includ- 
ing gravitational settling and radiative forces by Turcotte 
et al. (1998) (see Varenne & Monier 1999 for further de- 
tails). 

both field and open cluster evolved stars. First note that in 
classical models without rotation the variation in surface 
Li abundance is very similar for stars of different masses 
because it is only due to the first dredge-up dilution. More 
importantly, classical Li depletion starts at a too low T e ff 
(around ~ 5850 K) and leads to final Li abundances too 
high when compared to observational data. On the other 
hand, and even if no important lithium depletion is ex- 
pected at the surface of the main sequence stars on the hot 
side of the dip at the age of the Hyades, more lithium de- 
struction occurs inside rotating models compared to clas- 
sical ones. Indeed, due to the enlargement of the Li free 
radiative regions, important surface Li depletion occurs at 
much higher T c g than in classical models. Also, total Li 
depletion after dredge-up is more important in rotating 
models than in models without transport processes. The 
predicted evolution of the surface lithium abundance in 




Fig. 10. Lithium abundance versus effective temperature 
in open clusters (open symbols) and field (filled symbols) 
subgiants compared with our rotating models. Triangles 
indicate upper limits. Observations are from Lebre et 
al.(1999) (field), Wallerstein et al. (1994) (field), Pasquini 
et al. (2001) (NGC 3680), Gilroy (1989) (open clus- 
ters) and Burkhart & Coupry (1989,2000) (Hyades among 
other open clusters) . Dot-dashed lines represent results for 
standard models of 1.35M Q , 1.4M Q ,1.5M Q , 1.8M and 
2.2M Q . They all coalesce in a single line, standard dilution 
starting around 5750 K. Other lines are for models with 
rotation induced mixing and magnetic braking. Solid line, 
dotted line, long dashed line, dashed line and dot-dashed 
line are for 1.5M Q , 1.4M , 1.8M , 1.35M Q , 2.2M Q respec- 
tively. Refer to Tab. 1 for details concerning rotational 
velocities. 

our rotating models explains the behaviour and the disper- 
sion in field and open cluster evolved stars more massive 
than ~ 1.4 M Q . 

7. Conclusions 

To study the evolution of low mass rotating stars and the 
effects of rotation induced mixing on surface abundance 
patterns, we have presented here a first series of results 
concerning main sequence and subgiant Pop I stars lo- 
cated on the hot side of the lithium dip. We performed 
calculations including the treatment of transport of an- 
gular momentum and chemical species due to meridional 
circulation, shear turbulence and gravitational settling ac- 
cording to Maeder & Zahn (1998) formalism. In particu- 
lar, we take into account horizontal variations of mean 
molecular weight, and their action on the transport itself. 
These terms prove to be very important in defining the 
shape of the internal rotation profile. Indeed the horizon- 
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tal variations of (i and the derivative of Q, namely, A and 
8 respectively, have a mirrored behaviour, 6 following A 
when it is not null. As soon as the circulation is settled, the 
/z-currents grow larger and induce a stronger differential 
rotation. 

On the other hand, the effects of these terms on the 
surface abundances of light elements appear not to be con- 
straining compared with braking for instance. The major 
change introduced in the transport of chemicals when tak- 
ing /i-currents into account concerns the dominant con- 
tributor to the diffusion coefficient. In an "homogeneous" 
star (where there are no horizontal variations of /z), the 
effective diffusion coefficient D e ff is the major contribu- 
tor, while in the inhomogeneous case, turbulent diffusion 
coefficient dominates the transport of chemicals. 

Finally, contrary to what has been suggested in previ- 
ous work concerning rotation induced mixing in solid body 
rotating low mass stars (typically halo main sequence stars 
with M* < O.9M ), when considering extraction of angu- 
lar momentum through magnetic braking, the situation of 
11 creeping paralysis" of the circulation is not reached in the 
models presented here. It shall actually never be the case 
as long as differential rotation is allowed to exist for the 
mirrored terms never compensate each other. 

The point to be emphasized is the fact that these vari- 
ations of mean molecular weight are not to be neglected 
when considering rotating objects. They modify the action 
of mixing and can be of major importance after main se- 
quence evolution, for they constrain the internal rotation 
profile. 
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